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For every positive integer n let S(n) be the minimal positive integer m such 
that n | m!. For any positive number z > 1 let 


A(z) = = > S(n) a 


n<z 
be the average value of S on the interval {1, z]. In [6], the authors show that 
A(z) <c,z +2 (2) 


where cı can be made rather small provided that z is enough large (for example, 
one can take cy = .215 and co = 45.15 provided that z > 1470). It is interesting 
to mention that by using the method outlined in [6], one gets smaller and smaller 
values of cı for which (2) holds provided that z is large, but at the cost of increasing 
co! In the same paper, the authors ask whether it can be shown that 





2r 
A 
()< p= (3) 
and conjecture that, in fact, the stronger version 
T 
A(z) < EE (4) 


might hold (the authors of [6] claim that (4) has been tested by Ibstedt in the 
range z < 5. 10° in [4]. Although I have read [4] carefully, I found no trace of the 
aforementioned computation!). 


In this note, we show that is indeed the correct order of magnitude of 


A(z). 
For any positive real number z let (z) be the number of prime numbers less 
then or equal to z, 


x 
logz 


B(z)=zA(z)= X. Sín), (5) 
E(x) = 2.5 log log(z) + 6.2 + =. (6) 
We have the following result: 
Theorem. 
5(a(z) — r(vVz)) < A(z) < t(z)+ E(x) forall z>3. (7) 


Inequalities (7), combined with the prime number theorem, assert that 
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.5 < liminf —— = < oe ao) <A, 


= t z 











which says that is indeed the right order of magnitude of A(z). The natural 


z= 
log z 


conjecture is that, in fact, 
z T 
A(z) = — + o(—-). (8) 


Since 


T 1 3 
log z Tem —— > 
iogz (1+ Bogs) aay eo ee) for z > 59, 


it follows, by our theorem, that the upper bound on A(z) is indeed of the type (8). 
Unfortunately, we have not succeeded in finding a lower bound of the type (8) for 
A(z). 


The Proof 


We begin with the following observation: 


Lemma. 
Suppose that n = p?*...p?* is the decomposition of n in prime factors (we 
assume that the p;’s are distinct but not necessarily ordered). Then: 
l. 
S(n) < maxi- (2:P:)- (9) 
2. Assume that a,p, = max*_,(a;p;). [fai < pi, then S(n) = aipi. 


3. 
S(n) > a;(pi — 1) for alli=1, ..., k. (10) 


Proof. 

For every prime number p and positive integer k let ep(k) be the exponent at 
which p appears in k!. 

l. Let m > max*_,(a;pi). Then 


m m 
ep (m) = Be Z| >a; foris lh aeg k 
ps (m) 2 7 5 
This obviously implies n | m!, hence m > S(n). 
2. Assume that a; < pı. In this case, S(n) > œıpı. By 1 above, it follows that 
in fact S(n) = ayp1. 
3. Let m = S(n). The asserted inequality follows from 


ai < enlm) =) || < mye 


s>1 t s>1 
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The Proof of the Theorem. 


In what follows p denotes a prime. We assume z > 1. The idea behind the 
proof is to find good bounds on the expression 


B(z)- B(Vz)= > S(n). (11) 


Jrin<s 
Consider the following three subsets of the interval I = (Vz, z]: 


Cı = {n € I | S(n) is not a prime}, 

C2 = {n E€ I | S(n) =p < V3}, 

C3 = {n € I | S(n) = p > Vz}. 
Certainly, the three subsets above are, in general, not disjoint but their union covers 
I. Let 


D,(z) = > S(n) for i = 1, 2, 3. 
nec, 


Clearly, 
max(D;(z) | i= 1, 2, 3) < Bl) — B(VZ) < Di(z) + D(z) + Ds(2). (12) 


We now bound each D; separately. 


The bound for D;. 


Assume that m € Cı. By the Lemma, it follows that S(m) < ap for some 
p° || manda > 1. First of all, notice that S(m) < a,/m. Indeed, this follows from 
the fact that 

S(m) < ap<ap*??<a/m fora >2. 


In particular, from the above inequality it follows that p < ym < /z. Write now 
m = p%k. Since m < z, it follows that k < x/p*. These considerations show that 


Da Oe See Sar Py 
pyr 222 p pza? eyes) 


In the above formula (13), we used the fact that 


a= d 1 D 1 2 _ 2z— 2? 
ee t= —(-—)-1=(_) eee EES for |z| < 1 





with z = i/p. Since 





2p—1 5 
— for p > 3, 
(p — 1) ~ 4p üg 
it follows that 
5 Oo 1 1 
D(z) <2(3-24+2 D> -) =2(2.37541.25 JD =), (14) 
8 4 P P 
PEVE pSVZ 


2l 


From a formula from [5], we know that 


> = < loglog y + 1.27 for all y > 1. 
psy 


Hence, inequality (14) implies 
Dy(z) < 2(2.375 + 1.25 (loglog Vz + 1.27)) < 2(3.1 + 1.25 loglog z). (15) 


The bound for D» 
Assume that S(m) = p. Then m = py where p does not divide y. Since 
m > /z, it follows that 
Se 


—<y< 
p 


Since p < yz, it follows that at least one integer in the above interval isa multiple 
of p; hence, cannot be an acceptable value for y. This shows that there are at most 


BIR 





P| 
P 7 P 
possible values for y. Hence, 
z- yz 
Di(z)< > p-( ) < (£ - vz)n(v2). (16) 
p 
pSV 
Bounds for D3 
Assume S(m) = p for some p > yz. Then, m = py for some y < z/p. Hence, 
T 
Ds(z)= © p |5]. (17) 
Vi<psz 


Notice that, unlike in the previous cases, (17) is in fact an equality. Since z > |z] > 
.5z for all real numbers z > 1, it follows, from formula (17), that 


52(m(z) — x(V2)) < Da(z) < 2(x(z) — (V2). (18) 


Denote now by 
F(z) = 3.1+ 1.25 log log(z) 


From inequalities (12), (15), (16) and (17), it follows that 
52(4(z) — t(/z)) < Dg(z) < B(x) — B(/z) < Di(z) + Do(z) + D(z) < 


2F (2) + (z — V2)(Jz) + 2(m(z) — x(V2)) = 2x(2) — Ven(Vz) + F(z). (19) 


The left inequality (7) is now obvious since 


B(z) > B(Vz) + 5e(x(z) — #(Vz)) > 1 + .52(x(z) — (V2). 
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For the right inequality (7), let G(z) = zr(z). Formula (19) can be rewritten as 
B(z) - B(/z) < G(z) - (vz) + zF (2). (20) 


Applying inequality (20) with z replaced by yz, z'/*, ..., 24/2" until r? < 2 
and summing up all these inequalities one gets 


B(z) — B(1) < G(z) + y~ 21/2 Fy 1/2'), (21) 


s=0 


The function F(z) is obviously increasing. Hence, 


B(z) < 1+ G(z) + F(z) y oa (22) 


s=0 


To finish the argument, we show that 
EDD (23) 
i=l 


Proceed by induction on s. If s = Q, there is nothing to prove. If s = 1, this just 
says that z > /z which is obvious. Finally, if s > 2, it follows that z > 4. In 
particular, z > 2,/z or z — yz > yz. Rewriting inequality (23) as 


s 


b- & 


which is precisely inequality (23) for yz. This completes the induction step. Via 
inequality (23), inequality (22) implies 


B(z) <1+27(z) 4+ 22 F(z) = 1+ 27(z) +22 (3.1 +1.25loglogz) (24) 


or 
i 
A(z) < m(z) + m 6.2+2.5loglogz = a(x) + E(z). 


Applications 


From the theorem, it follows easily that for every € > 0 there exists zgo such 


that 
r 


A(z) < (14+ eee 


(25) 





In practice, finding a lower bound on Zo for a given €, one simply uses the theorem 
and the estimate 


z ff 3 
m(z) < EER (1 + Toga) for z > 1. (26) 
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(see (5]). By (7) and (26), it now follows that (25) is satisfied provided that 


is (Š + 5(2)). 


loge ~ €\2iog?z 


For example, when € = 1, one gets 


ž 
A(z) < E for z > 64, (27) 
for € = .5, one gets i 
A(z) < lS for z > 254 (28) 
and for € = 0.1 one gets 
A(z) <1.1—— forz > 3298109. (29) 
logz 


Of course, inequalites (27)-(29) may hold even below the smallest values shown 
above but this needs to be checked computationally. 
In the same spirit, by using the theorem and the estimation 


z 1 
R m > 
a(z) > eee (1 + ee) for z > 59 


(see [5]) one can compute, for any given €, an initial value zo such that 


T 








A(z) > (.5- Nees for z > Zo. 
For example, when ¢ = 1/6 one gets 
A(z) > ae for z > 59 (30) 
3 logz ree 


Inequality (30) above is better than the inequality appearing on page 62 in [2] which 
asserts that for every a > 0 there exists zg such that 


A(z) > z% for z > Zo (31) 
because the right side of (31) is bounded and the right side of (30) isn’t! 
A diophantine equation 


In this section we present an application to a diophantine equation. The ap- 
plication is not of the theorem per se, but rather of the counting method used to 
prove the theorem. 

Since S is defined in terms of factorials, it seems natural to ask how often the 
product S(1) - S(2) -...- S(n) happens to be a factorial. 

Proposition. 


The only solutions of 


S(1) - S(2)-...-S(n) = m! (32) 
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are given by n= m € {1, 2, ..., 5}. 


Proof. 

We show that the given equation has no solutions for n > 50. Assume that this 
is not so. Let P be the largest prime number smaller than n. By Tchebysheff’s 
theorem, we know that P > n/2. Since S(P) = P, it follows that P | m!. In 
particular, P < m. Hence, m > n/2. 

We now compute an upper bound for the order of 2 in S(1) - S(2)-...- S(n). 
Fix some J > 1 and assume that k is such that 2° |j S(k). Since 


S({k) = max(S(p*) | p* || k), 


it follows that 2° || S(p%) for some p° || k. 
We distinguish two situations: 


Case 1. 

p is odd. In this case, 2fp | S(p*). If 8 = 1, then a = 2. If 8 = 2, then a = 4. 
For 6 > 3, one can easily check that a > 29 — 8 + 1 (indeed, if a < 2° — 2, then one 
can check that p° | (2°p — 1)! which contradicts the definition of S). In particular, 
p? -2+ | k. Since 27-1 > x +1 for x > 3, it follows that a > 22-1 +2. Since 
k < n, the above arguments show that there are at most 


a for f =1, 2 
and 
n 
p2 for B >3 


integers k in the interval [1, n] for which p | k, S(k) = S(p%), where a is such that 
p° || k and 2° || S(k). 
Case 2. 
p= 2. If J = 1, then k = 2. If 8 = 2, then k = 4. Assume now that 8 > 3. 
By an argument similar to the one employed at Case 1, one gets in this case that 
a > 28 — B. Since 2° || k, it follows that 92°-8 | k. Since k < n, it follows that 


there are at most N 


227- 
such k’s. 


From the above anaysis, it follows that the order at which 2 divides S(1)-S(2)- 
.- S(n) is at most 


a<d+n D (3 k q? pe eq) + Las (38) 


PAR 
p odd 


(the number 3 in the above formula counts the contributions of S(2) = 2 and 
S(4) = 4). We now bound each one of the two sums above. 
For fixed p, one has 


1 2 3 4 y p? 
+44 Spigtutatmt-<La Gop 
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Hence, 


1 2 p p 
~+—+) =! < oy < -245 40) 
Z (z p’ 2 p 7) p (r 7 1 


We now bound the second sum: 


B _3, 4 5 3 B _ 
De sep = ostotgrt.-< gt 5 moa 
p>3 B>3 
¥+2 . 1/15 31 
+ 167 p) = 56 +7 (i + as) < 09 (41) 
From inequalities (38), (40) and (41), it follows that 
e2 < 34+ .344n. (42) 


We now compute a lower bound for ez. Since e2 = e2(m!), it follows, from Lemme 
1 in [1] and from the fact that m > n/2, that 


log(m+1) _ n  log(n/2+1) 
> eee —EE > — — ba ete 
2am log2 ^? log 2 (43) 


From inequalities (42) and (43), it follows that 


ETTE E o 
= log 2 
which gives n < 50. One can now compute S(1) - S(2)-...- S(n) for all n < 50 


to conclude that the only instances when these produci? are factorials are n = 


i E ee, 
We conclude suggesting the following problem: 


Problem. 
Find all positive integers n such that S(1), S(2), ..., S(n*) can be arranged in 
a latın square. 


The above problem appeared as Problem 24 in SNJ 9, (1994) but the range 
of solutions was restricted to {2, 3, 4, 5, 7, 8, 10}. The published solution was 
based on the simple observation that the sum of all entries in an n x n latin square 
has to be a multiple of n. By computing the sums B(x?) for z in the above range, 
one concluded that B(z?) Æ 0 (mod z) which meant that there is no solution for 
such z’ses. It is unlikely that this argument can be extended to cover the general 
case. One should notice that from our theorem, it follows that if a solution exists 
for some n > 1, then the size of the common sums of all entries belonging to the 
same row (or column) is © na(n?). 


Addendum 


After this paper was written, it was pointed out to us by an annonymous referee 
that Finch [3] proved recently a much stronger statement, namely that 


lim reste) A(z yas - = 0.82246708... (44) 


T=% 
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Finch’s result 1s better than our result which only shows that the limsup of the 
expression log(z)A(z)/z when z goes to infinity is in the interval [0.5, 1]. 
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